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Abstract 


The  evolution  of  unsteady  boundary  layers  on  oscillating  airfoils  Is 
Investigated  by  solving  the  governing  equations  by  the  Characteristic  Box 
scheme.  The  difficulties  associated  with  computing  the  first  profile  on  a 
given  time-line,  and  the  velocity  profiles  with  partial  flow  reversal  are 
solved.  A  sample  calculation  has  been  performed  for  an  external  velocity 
distribution  typical  of  those  found  near  the  leading  edge  of  thin 
airfoils.  The  results  demonstrate  the  viability  of  the  calculation 
procedure. 


Introduction 


The  effect  of  unsteady  motion  of  an  airfoil  on  Its  stall  behavior  Is  a 
problem  of  great  importance  In  many  types  of  fluid  motion,  for  example, 
helicopter  rotor  blades  and  jet  engine  compressors.  A  recent  detailed 

study  has  been  published  by  Carr,  McAlister  and  McCroskey^ ;  this  study 
reveals  a  very  complicated  phenomenon  which  depends  In  a  subtle  way  on  a 
large  number  of  parameters.  There  Is  one  Important  characteristic  of  the 
data  they  have  compiled  which  serves  as  a  focal  point  of  the  approach 
which  we  are  planning  to  apply  sometime  later  In  our  studies;  namely  that 
at  some  stage  of  the  cycle,  a  large  vortex  Is  formed  near  the  surface  of 
the  airfoil  and  very  shortly  afterwards  stall  occurs.  It  seems  also  that 
the  stall  Is  associated  with  flow  reversals  In  the  unsteady  boundary  layer 
which  may  spread  downstream  or  upstream  depending  on  the  leading-edge 
radius  of  the  airfoil. 

The  present  paper  concerns  Itself  with  the  calculation  of  boundary-layer 
characteristics  of  an  oscillating  airfoil  In  order  to  Investigate  the 
evolution  of  unsteady  boundary  layers  on  such  airfoils.  It  Is  one  phase 
of  a  study  which  will  be  extended  later  In  the  hooe  of  throwing  light  on 
the  dynamic  stall  problem.  In  our  present  study  we  focus  our  attention  on 
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the  calculation  of  time-dependent  boundary  layers  for  a  given  pressure 
dl strlbutlon. 

There  are  three  difficulties  associated  with  unsteady  boundary  layers  that 
requires  careful  attention.  First  of  all  appropriate  initial  values  at 
t>0  rust  be  chosen  for  the  velocity  distribution.  Strictly  speaking  they 
can  be  arbitrary  but  in  that  event,  the  values  of  3u/3t  at  t  ■  0  Is 
non-zero  and  this  Implies  an  Invlscld  acceleration  of  the  fluid  In  the 
boundary  layer  and  In  consequence  a  velocity  of  slip  begins  to  grow  at  the 
wall.  This  is  smoothed  out  by  an  Inner  boundary  layer  Initially  of 
thickness  (vt)*^  In  which  viscous  forces  are  of  Importance.  Thus  a 
double  structure  develops  In  the  boundary  layer  which  may  be  treated  by  a 
generalization  of  the  Keller-box  scheme^.  However,  If  interest  Is 
centered  on  the  solution  at  large  times,  this  feature  may  be  reduced  In 
importance  by  requiring  that  the  Initial  velocity  distribution  satisfies 
the  steady-state  equation  with  the  instantaneous  external  velocity.  In 
addition  it  is  necessary  to  smooth  out  the  external  velocity  ug(x,t)  so 
that  aue/dt  >0  at  t  *  0  and  then  standard  methods  may  be  used  and 
are  stable.  The  use  of  the  smoothing  function  makes  for  some  loss  of 
accuracy  at  small  values  of  t  but  the  error  soon  decays  to  zero  once  the 
required  value  of  ug  is  specified.  In  the  present  problem  the  choice 
of  parameters  in  the  specified  external  velocity  distribution  is  such  that 
the  smoothing  function  Is  actually  unnecessary  but  it  can  easily  be 
Incorporated  Into  the  scheme. 

The  second  difficulty  arises  when  u  changes  sign  over  part  of  the 
profile  at  some  x-statlon  where  the  x-axIs  Is  parallel  to  the  streanvlse 
direction  and  u  Is  the  corresponding  velocity  conponent.  Normally  this 
does  not  occur  and  one  can  Integrate  away  from  the  profile  In  the 
direction  of  positive  u  without  any  difficulty  by  using  a  standard 
numerical  scheme.  However,  If  the  change  In  sign  of  u  does  occur,  we 
encounter  numerical  Instabilities  since  In  the  negative  u-reglon  we  would 
be  Integrating  against  the  stream.  The  Instability  can  be  avoided  by 
changing  the  scheme  either  to  the  zig-zag  box  or  the  characteristics  box. 
These  new  schemes  have  already  been  shown  to  be  effective  In  such 
circumstances  when  the  flow  Is  unsteady^  and  In  three-dimensional 
flows<4-6>.  The  essence  of  these  schemes  Is  that,  to  an  Increasing 
extent,  they  take  Into  account  the  fact  that  small  disturbances  are 
carried  along  with  the  local  fluid  velocity. 
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The  third  difficulty  arises  when  It  Is  desired  to  compute  the  first 
velocity  profile  at  the  new  time  line.  Given,  as  we  are,  the  complete 
velocity  profile  distribution  on  the  previous  time  line,  there  Is  In 
principle  no  difficulty  In  computing  values  on  the  next  time-line  by  an 
explicit  method,  but  If  we  wish  to  avoid  the  stability  problems  associated 
with  such  a  method  by  using  an  Implicit  method,  we  are  Irmedlately  faced 
with  the  problem  of  generating  a  starting  profile  on  the  new  time-line. 


In  order  to  explain  the  problem  further.  It  Is  Instructive  to  see  what 
happens  to  the  stagnation  point  as  a  function  of  time.  For  this  purpose 
let  us  assume  that  the  external  velocity  distribution  for  an  oscillating 
airfoil  Is  given  by  the  following  function. 


u 


x  ♦  £0  (1  +  A  sin  «t) 


e 


(1) 


(1  +  x2),s 

where  A  and  zQ  denote  parameters  that  need  to  be  specified.  This 
equation  is  a  good  approximation  to  the  external  velocity  distribution 
near  the  leading  edge  of  a  class  of  thin  airfoils  at  variable  angles  of 
attack  and,  when  A  »  0,  has  recently  been  used  by  Cebecl,  Stewartson,  and 

Williams^  to  study  leading-edge  separation  In  steady  flow. 


Since  by  definition 
Is  given  by 


u  *  0  at  the  stagnation  point.  Its  location,  x  , 

0  5 


X  »  (1  ♦  A  sin  ut) 


(2) 


and  so  the  upper  and  lower  surfaces  of  the  airfoil  as  functions  of  time 
are  defined  In  particular  by 
us  take  A  •  i, 

shown  In  Figure  1  for  one  cycle  (0  <  t  ^  8).  When  t  •  2,  the  stagnation 
point  xs  Is  at  -2  C0,  when  t  *  6  It  Is  at  0,  etc.  If 


x  >  x  and  x  <  x  .  For  example,  let 

5  5 

u  «  ir/4  and  plot  x$/£q  In  the  (t,x)  plane,  as 


x$  were 


fixed  we  could  assume  that  u  ■  0  at  x  »  x$  for  all  time  and  all 
Further  profiles  at  this  time-line  then  follow  by  use  of  one  of  the  box 
schemes  that  have  been  developed.  However  x  Is  not  fixed  and  It  Is 
clearly  unjustified  to  assume  a  priori  that  u  ■  0  there 
therefore  we  use  the  characteristic  box,  with  an 
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K 

JL 


extrapolated  normal  velocity,  to 
complete  the  profile  at  x  ■  Xj,  the 


nearest  x  station  to  x$  on  the 
new  time-line.  Thus  we  avoid  the  need 


-2.0  -1.0  0  1, 
xs/Co 

Fig.  1.  Variation  of  stagnation 
point  with  time  for  one  cycle 
according  to  Eq.  (2),  with 
w  *  ir/4. 


Governing  Equations 


to  use  two  x-statlons  at  the  new  time¬ 
line  to  dlmenslonallze  the  governing 
equations.  Once  this  step  Is 
completed,  the  determination  of 
further  profiles  on  the  new  time  line 
proceeds  normally  and  moreover  we  can 
Improve  our  Initial  estimate  of  the 
normal  velocity  at  x  ■  x^  by 
Iteration. 


The  governing  boundary-layer  equations  for  an  incompressible  laminar  or 
turbulent  flow  past  an  oscillating  airfoil  are  well  known  and,  with  the 
eddy  viscosity  (em)  concept,  they  can  be  written  as 

2H.  +  2I.0  (• 

ax  ay  v 


M  +  u  iH.  +  v  — 

at  ax  ay 


It  +ue 


Jx  ay 


where  t  *  v(3u/3y)  -  It.  In  the  absence  of  mass  transfer,  these 
equations  are  subject  to  the  boundary  conditions  given  by 

y-0,  u  ■  v  ■  0;  y  «,  u  ♦  ue(x,t)  (5) 

The  presence  of  the  Reynolds  shear  stress  term,  -  u' v' ,  requires  a 
closure  assumption;  In  our  study  we  use  the  algebraic  eddy -viscosity 
formulation  developed  by  Cebecl  and  Smith.  For  details,  see  ref.  8. 

To  complete  the  fornulatlon  of  the  problem.  Initial  conditions  must  be 
specified  In  the  (t,y)  plane  at  some  x  •  x0  either  on  the  lower  or 
upper  surface  of  the  airfoil  (see  Fig.  I)  as  well  as  Initial  conditions  In 
the  (x,y)  plane  on  both  surfaces  of  the  airfoil.  In  the  latter  case.  If 
we  assume  that  steady-flow  conditions  prevail  at  t  »  0,  then  the  initial 
conditions  In  the  (x,y)-p1ane  can  easily  be  generated  for  both  surfaces  by 


-3 ■■■  :%m 
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solving  the  governing  equations  for  steady  flow,  which  In  this  case,  are 
given  by  Eq.  (3)  and  by 


u  i!L  +  v  iH.  ■  u  ^£  +  J_(bJJi) 

ax  3y  e  dx  ay  v  ay  ' 


(6) 


where  b  »  v  +  cm.  There  Is  no  problem  with  the  initial  conditions  for 

R1 

Eqs.  (1)  and  (6)  since  the  calculations  start  at  the  stagnation  point. 


Generation  of  the  initial  conditions  in  the  (t,y)  plane  at  x  »  xQ, 
which  Is  one  of«the  purposes  of  our  study*  is  not  so  easy  as  was  discussed 
In  the  previous  section.  The  following  section  describes  our  solution 
procedure  at  t  »  0  and  t  >  0. 


Solution  Procedure 


Me  use  Keller's  Box  method  to  solve  the  governing  equations  of  the 

previous  section.  This  is  a  two-point  finite-difference  method  which  has 

been  used  to  solve  a  wide  range  of  parabolic  partial -differential 

g 

equations  as  discussed  by  Bradshaw,  Cebecl  and  Whltelaw  .  The  solution 
procedure  for  t  *  0  (Eqs.  (3)  and  (6)) and  t  >  0  (Eqs.  (3)  and  (4))  are 
described  separately  below. 


Solution  Procedure  for  t  ■  0 

As  explained  In  the  Introduction  any  velocity  profile  may  be  chosen  at 
t  *  0  to  Initiate  the  computation  but  It  is  convenient  to  select  one 
which  obviates  the  need  for  double-structured  numerical  schemes  and  joins 
smoothly  on  to  the  solution  for  t  >  0.  Me  Insist  that  these  profiles 
satisfy  the  steady  state  equations  with  ue  given  by  (1)  and  t  ■  0. 

The  details  of  the  procedure  for  computing  these  profiles  differs  slightly 
from  previous  procedures  used  In  steady  two-dimensional  flows  where  we 
have  used  the  definition  of  the  stream  function  and  reduced  Eqs.  (3)  and 
(6)  to  two  first-order  ordinary  differential  equations  and  to  one 
partial -differential  equation.  Here  we  consider  the  solution  of  Eqs.  (3) 
and  (6)  without  the  use  of  the  stream  function.  For  this  purpose,  with 
primes  denoting  differentiation  with  respect  to  y,  we  let 
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u*  -  f 


(7a) 


and  write  Eqs.  (3)  and  (6)  as 


■  3u 
v  »  -  — 

ax 


(7b) 


(bf)'  -  fv 


d 

3x 


(iyuf) 


(7c) 


The  finite-difference  approximations  of  Eqs.  (7)  are  also  somewhat 
different  than  those  reported  In  our  previous  studies  dealing  with 
two-dimensional  flows^.  All  quantities  except  for  the  normal  velocity 
component  v,  are  centered  at  the  center  of  the  box  (y 


J - 1/2  * 

x.  1/2^’  see  ^9*  by  tak-in9  the  values  of  each  parameter,  say 
the  four  corners  of  the  Box,  that  Is, 


q,  at 


1-1/2 


•r«! 


,1-1/2, 


i 

4j-i 


+  q<  ,  +  q< 


i-i; 

4j-r 


(8a) 


j -1/2 

However,  the  centering  of  the  y-veloclty  component  v  Is  done  by  writing 
it~as 


1-1/2 


1 

7 


(v]-i/2  ♦  v]:}/2) 


(8b) 


1-1/2  „ 
Vj  so 


- 4 - 9 - f  "j-l 


l 


T 


j-l/2 

The  unknown  parameters  In  Eqs.  (8)  correspond  to  qj  and 
that  when  a  solution  of  the  system  given  by  Eqs.  (7)  Is  obtained,  f  and 
u  are  computed  at  (1,j)  and  v  at  (1-1/2,  j).  This  modified 

centeHnq  procedure  Is  necessary  in 
order  to  avoid  oscillations  due  to  the 
use  of  the  continuity  equation  In  the 
form  given  by  Eq.  (7b)  rather  than  the 
use  of  the  stream  function.  The 
centering  of  Eqs.  (7)  and  the 
subsequent  linearization  procedure  by 
Newton's  method  allows  the  resulting 
linear  system  to  be  written  In  the 
form<10> 


y(j) 

yj 

yj-*i 

yj-i 


J 


M-l  *i-h 


ll 


x(  1 ) 


Fig.  2.  Net  rectangle  for 
finite-difference  approximations. 


*J  ’  6uj-l  *  -¥•  (5fj  +  ifj-l]  “  (r3}j-l 


(9a) 


6- 


b^t  U’j  -  *Vi>  *  <*7>j  ({Uj  *  {Uj-i>  ■  <ri>j 

‘sl>jSfJ  *  <s2>jJ,J-l  *  <s3>JS,J  *  <$4>j{,j-l  *  <sS>35u3 

+  (sg)j4Uj-l  "  (r2}j 

After  linearization,  the  wall  boundary  conditions 


(9b) 


(9c) 


u  »  0,  v  »  0  (10a) 

0  0 

and  the  edge  boundary  condition 

u.  *  u  (10b) 

J  6 

become 

5u .  *  0,  5V  *  0,  6u .  «  0  (11) 

O  O  J 


The  resulting  linear  system  consisting  of  those  given  in  Eqs.  (9)  and  by 
those  in  Eg.  (11)  can  be  solved  by  the  block  elimination  method  discussed 
by  Cebeci  and  Bradshaw' 


Solution  procedure  for  t  >  0 


We  have  already  pointed  out  that  if  au  /at  /  0  when  t  ■  0,  a 
double  structure  scheme  should  strictly  be  used  to  advance  the  solution 
from  t  »  0.  However  the  choice  of  parameters  In  our  study  Is  such  that 
aue/at  Is  small  and  the  difficulties  that  arise  from  using  a  standard 
method  are  of  a  sufficiently  minor  nature  that  no  further  refinement  Is 
necessary.  For  larger  values  of  the  relevant  parameters  It  Is  easy  to 

Incorporate  a  smoothing  function  Into  u„  and  one  can  always  use  the 

(21  e 
general  method' 


Nevertheless  there  Is  still  the  difficulty  about  obtaining  the  velocity 
profile  on  the  first  x-statlon  at  any  new  time-line.  It  can  be  resolved 
with  the  use  of  the  characteristic  box  method  developed  by  Cebeci  and 
Stewartson^).  Defining  the  streamline  by 


(12) 
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and  using  the  definition  of  f',  and  with  s  denoting  the  local 
streamline,  we  write  Eq.  (4)  as  _ 

t t(t|  |  |  '  |  ,  ' ,  <»fi  -  fv  W1  *  u«  is* 

ViT*“™ u«  .  JTT7  03) 

\  To  obtain  the  solution  of  the  unsteady 

k  - j - boundary. layer  equations  given  by  Eqs. 

v/  p  JpE i  \  at  (7a, b)  and  (13)  at  the  first  x-statlon 

k*1  - - * - * — * —  — on  either  side  of  the  stagnation  line, 

1+1  let  us  consider  the  grid  of  Fig,  1  and 

difference  molecule  for  the  direct  our  attention  to  the  point 

Characteristic  Box  2.  denoted  by  1  (see  Fig.  3). 


To  write  the  difference  approximations  of  Eq.  (13)  we  define 


as.  -  At/cos  0i  (14) 

J  J 

where,  with  u.  denoting  an  average  velocity,  we  compute  a,  from 
J  J 

a .  *  tan  ^  u  .  ( 15) 

J  J 

assuming  that  at  first  v  at  point  P  Is  known  and  Is  equal  to  Its  value 
i-1/2  k-1 

at  v  .  This  assumption  decouples  the  continuity  equation,  Eq. 

(7b)  from  Eqs.  (7a)  and  (13)  and  reduces  the  problem  to  a 
“two-dimensional"  one  with  f  and  u  being  the  only  unknowns.  The 
finite-difference  approximations  of  Eq.  (7a)  are  written  In  the  usual  way 
and  the  finite-difference  approximations  to  Eq.  (13)  are  written  by 
centering  It  at  point  P.  This  procedure  leads  to 


. -1  , .1 ,k  ,1  ,k.  1  ,k  . 

Vi(v  •  yp  ■  “j-i/2  ■  ° 


(16a) 


(bv)i*k-(bv)]:;  (bv^^-^v);:;-1  .  P 

- - 7  Lf j-1/2  +  rj-l/2JLvj-l/2J 

.  .  -  ,  1  ,k  n,k-l 

,  i - r1,k  i - *n,k-l  u  *  -  u  ’ 

+  \  [(VlTuf)  +  Ofl  ♦  U2)  ]  (  e-  ~  - ) 

J 

mird77)itk  *  o/r77)n,k'13 

2  1  J-1/2  ,j-l/2J  As.  . 


j  (16b) 
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The  profiles  bj,  fj  and  u^  as  well  as  u^  at  (n,k-l)  are 

obtained  by  Interpolating  the  profiles  at  (1,k-l)  and  (1-1, k-1).  To 

find  the  angle  «. ,  we  define  u.  In  Eq.  (15)  by 
J  J 

<17> 

Since  the  system  given  by  Eqs.  (16)  Is  linear,  there  Is  no  need  for 
linearization  and  we  solve  it  subject  to  the  two  boundary  conditions, 
namely, 

uo  *  °»  uj  *  ue  (18) 

by  using  the  block-elimination  method  In  which  case  the  matrices  are 
2x2.  We  shall  refer  to  this  scheme  as  Characteristic  Box  2. 

Once  a  solution  of  Eqs.  (7a)  and  (13)  has  been  obtained,  we  compute  v 
from  Eq.  (7b)  which,  in  finite-difference  form,  for  the  center  of  the  net 
rectangle,  point  E,  can  be  written  as  (see  Fig.  3) 


Here  Vj  denotes  the  value  of  v^  at  E  and  u^  is  given  by  Eq. 

(17).  Since  the  right-hand  side  of  Eq.  (19)  Is  known,  we  can  solve  this 

equation  for  V.  and  with  V  *  0,  find  V.  for  1  <  j  <  J.  We 
J  o  J  —  p 

then  substitute  this  new  value  of  into  Eq.  (16)  for  v  and 

solve  the  system  again  to  compute  new  values  of  Vj.  This  procedure  Is 
repeated  until  convergence. 

For  convenience  we  use  the  same  procedure  to  compute  point  2  to  the  left 

of  point  1.  Once  two  points  on  a  given  t-llne  are  computed  by  this 

procedure,  we  then  use  the  values  of  V<  at  E.  and  Et,  compute  a 
p  J  z  l 

new  value  v  and  repeat  the  solution  procedure  for  Eqs.  (7a)  and  (13), 
and  later  Eq.  (19).  After  that  the  stations  to  the  left  of  point  2  and 
the  stations  to  the  right  of  point  1  are  computed  by  using  the  Regular  Box 
scheme  If  there  Is  no  flow  reversal  across  the  layer  and  by  the  Character¬ 
istic  Box  scheme  If  there  Is  flow  reversal.  The  "new"  Characteristic  Box 
scheme  Is  now  slightly  different  than  the  Characteristic  Box  2  so  we  shall 
refer  to  It  as  Characteristic  Box  3. 


To  describe  the  Characteristic  Box  3  scheme  which  solves  Eqs.  (7a, b)  and 
(13)  without  decoupling  the  continuity  equation  from  Eqs.  (7a)  and  (13), 
we  consider  the  sketch  shown  in  Figure  4. 
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Fig.  4.  Notation  for  Characteristic  Box  3. 


Using  the  Zig-Zag  Box  scheme  discussed  In  detail  in  Ref.  5,  we  write  Eq. 
(7b)  in  the  following  finite-difference  form 


"jii  «»3  -  *i-i»  * 


where 


®  "  Tx 
h  1 


(x1+l  - 


1 


1+1  -  x1-l>  (xi  -  x1-l>  ’ 

1  /..I  »k-l  1+1, k-l 


xi  •  xi-i 


x1+l  -  XM 


7 


(uj-l/2  "  u j - 1/2 


) 


Since 


’1-1/2  V1 -3/2  Txi_3^2  -  xi .5/2 

The  relation  between  ^  and  ^  can  written  as 


» (xm/2  ~  Ibllil  rv  _v  ) 

^ r  '  1-3/2  vi-5/2' * 


vi  *  ^  +  V  Vl/2  +  6  3 

where  b2  and  B3  are  9lven  by 


(20) 


(21) 


(22) 


(xj  “  x1-l/2^ 

'  '  ‘“1-3/2  -  *t-l/2>  ’  *3  * 

Introducing  Eq.  (22)  Into  Eq.  (20)  and  rearranging,  we  get 

hj-l  ‘  »2><5J  -  5J.1>  *  ^uj-l/2  ■  «4 

where 

64  -  ♦  ®ujli/2k  +  hj-l  1  "  ®3^ 


(23) 


(24) 

(25) 
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As  In  Characteristic  Box  2,  we  center  Eq.  (13)  at  the  midpoint  of  (1,k) 
and  (1»n)  to  get  the  finite-difference  equations  given  by  Eq*  (16b)  with 
vP  being  obtained  by  linear  Interpolation  of  and 

which  Is 


VP  •  Vj  +  (Xc  -  Xj)  - - 


xn  -  xl 


(26) 


Equations  (16a),  (24)  and  (16b)  are  then  linearized  by  Newton's  method, 
and  again  are  solved  by  the  block  elimination  method. 


When  there  Is  no  flow  reversal  across  the  layer,  we  use  the  Regular  Box 
scheme  described  In  detail  In  Ref.  5. 


Results  and  Discussion 

To  date  calculations  have  been  carried  out  in  only  one  test-case,  namely 
when  g  ■  0.10,  A  •  1,  «  ■  »/4,  and  for  a  limited  range  of  x 
(  | x |  <  0.3).  With  the  use  of  the  various  procedures  described  in  this 
paper  the  calculations  were  quite  straightforward  and  the  formal  validity 
and  efficacy  of  the  numerical  schemes  were  established.  The  results  are 
summarized  In  Figs.  5-7.  In  Fig.  5  we  display  the  variation  of  wall-shear 
with  time  at  different  x-statlons  and  In  Fig.  6  the  variation  with  x  at 
different  times.  These  graphs  are  entirely  In  line  with  expectations  and 
we  note  that  the  flow  reversal  at  the  wall  Is  quite  smooth.  A  similar 
remark  applies  to  the  velocity  profiles  on  either  side  of  the  stagnation 
line  displayed  In  Fig.  7. 

The  next  phase  In  our  studies  Is  to  extend  the  computations  to  larger 
values  of  A,£0  and  smaller  values  of  <■>  so  as  to  more  closely 
approach  the  conditions  of  dynamic  stall.  It  Is  of  Interest  to  comment  on 
the  fluid  mechanical  problems  that  may  then  arise.  First,  If  £0  Is 
Increased  beyond  1.155,  the  steady-state  solution  at  t  ■  0  separates  on 
the  upper  side  of  the  airfoil  and  the  calculation  terminates.  This  Is  not 
a  serious  drawback  unless  the  unsteady  boundary  develops  a  singularity 
because  the  smoothing  function  mentioned  earlier  may  be  adapted  to  ensure 
that  £Q  is  Initially  less  than  1.155  and  rises  to  a  value  greater  than 
that  after  a  finite  time.  The  unsteady  boundary  layer  then  Includes 
regions  of  reversed  flow  which  may  well  become  extensive  If  A  Is  also 


allotted  to  Increase  to  mimic  more  closely  the  conditions  of  dynamic 

stall.  Even  If  the  boundary  layer  remains  smooth,  the  displacement 

thickness  may  then  become  much  thicker  and  have  a  significant  modifying 

effect  on  the  external  flow.  It  would  be  useful  then  to  consider  an 

Interactive  problem  In  which  the  external  stream  depends  In  part  on  the 

displacement  thickness  thus  generalizing  the  studies  reported  for  steady 

flow*7*.  Consideration  then  has  to  be  given  to  the  variation  of 

circulation  with  time  which  may  lead  to  a  more  complicated  expression  for 

the  dependence  of  u  on  the  displacement  thickness  than  was  used* In 

6 

Ref.  7  but  the  computation  should  not  be  any  more  complicated  as  a 

result.  Finally,  In  order  to  mimic  the  dynamic  stall  problem  most 

effectively^*’**^ ,  <■>  should  be  reduced  to  very  small  values  (as 

typical  of  dynamic  stall  problem).  So  long  as  «  >  0,  the  difficulties 

reported  in  Ref.  3  at  separation  In  uninteracted  flows  and  In 

post-separation  flows  otherwise,  should  not  be  present.  On  the  other 

(121 

hand',  van  Oommelen  and  Shen'  '  have 
provided  quite  strong  evidence  that  a 
singularity  can  occur  In  an  unsteady 
boundary  layer  for  which  the  external 
velocity  Is  steady.  This  phenomenon 
Is  still  somewhat  controversial ^*^ 
but  there  seems  no  doubt  that  the 
boundary  layer  will  exhibit  dramatic 
properties  for  small  enough  values  of 
-0.2  -o.i  o  o  o  *°-1  0  °-1  °-2  °’J  «  and  It  Is  possible  that  these  may 

Fig.  5.  Variation  of  wall  shear  »,,e  ,urth*r  ,ns1»ht  ,nt0 
parameters  with  time  at  differ-  stall, 
ent  x-statlons. 


Fig.  6.  Variation  of  skin-friction  coefficient  with  x  at  different 
t-lntervals.  Note  t  ■  8  corresponds  to  one  cycle. 
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Fig.  7.  Velocity  profiles  In  the  Immediate  neighborhood  of  the  stagnation 
lines  at  different  times,  u  *  0  on  the  dashed  line  at  the  specified  time 
and  serves  to  "define"  upper  or  lower  surfaces. 
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